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Abstract 

We analyse the effect of quenched uncorrelated randomness cou- 
pling to the local energy density of a model consisting of N coupled 
two-dimensional Ising models. For N > 2 the pure model exhibits 
a fluctuation-driven first order transition, characterised by runaway 
renormalisation group behaviour. We show that the addition of weak 
randomness acts to stabilise these flows, in such a way that the tra- 
jectories ultimately flow back towards the pure decoupled Ising fixed 
point, with the usual critical exponents a = 0, u = 1, apart from 
logarithmic corrections. We also show by examples that, in higher 
dimensions, such transitions may either become continuous or remain 
first order in the presence of randomness. 
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The effect of quenched randomness coupling to the local energy density 
of a system which, in its absence, undergoes a continuous phase transition is 
well understood from the point of view of the renormalisation group version 
of the Harris criterion |lj. When the specific heat exponent a of the pure 
model is negative, weak randomness is irrelevant from the renormalisation 
group point of view, and the pure fixed point is stable. On the other hand, 
when a > it is relevant, and, at least when the cross-over exponent a is 
small, it may be argued that the critical behaviour is controlled by a new, 
random, fixed point close by. 

The effect on systems which undergo thermal first order transitions is 
more dramatic. It was argued some time ago by Imry and Wortis[[| that, in 
two dimensions, such systems should always exhibit a continuous transition 
in the presence of such randomness. This is because the random impurities 
couple to the local energy density in much the same way that a random 
field couples to the local magnetisation in an Ising system. In dimensions 
d < 2, the Imry-Ma argument |§ implies that such random fields should 
destroy the ordered phases at low temperature, and therefore also the first 
order phase boundary between them. A similar argument, applied to ran- 
domness coupling to the local energy density, then implies that a non-zero 
latent heat is impossible in two dimensions in random systems whose pure 
versions exhibit such behaviour. This argument has been rediscovered and 
put on a rigorous basis by Aizenman and Wehr[|J, and is supported by the 
phenomenological and approximate renormalisation group arguments of Hui 
and Berker||. Monte Carlo work of Chen, Ferrenberg and Landau || on the 
q = 8 state Potts model and of Domany and Wiseman on the Ashkin- 
Teller and 4-state Potts models supports this conclusion, and goes further: 
the continuous transition found by these workers exhibits critical exponents 
which are consistent with those of the pure Ising model, namely j/u ~ 1.75, 
(3/v ~ 0.125 and a ~ 0. An argument explaining these findings has been put 
forward by Kardar, Stella, Sartoni and Derridag]. They study the properties 
of an interface in the g-state random-bond Potts model at low temperatures. 
For q ^ 2 this has a branching structure, but the authors argue, on the basis 
of simplified recursion relations which are exact on a hierarchical lattice, that 
the critical behaviour where the interfacial free energy vanishes is governed 
by a zero-temperature fixed point (as in the random field problem) , and the 
Widom exponent \i which governs the vanishing of the surface tension is in- 
dependent of q for sufficiently large q, being numerically consistent with the 
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Ising value fi = 1. 

The first order transition in the pure g-state Potts model is of mean field 
type, that is it is already predicted on the basis of mean field theory. Such 
'strong' first order transitions are described within the renormalisation group 
by zero-temperature, discontinuity fixed points, characterised by a relevant 
renormalisation group eigenvalue y = d whose scaling field couples to the 
local energy density Quenched randomness coupling to this has eigenvalue 
d — 2(d — y) = d, and is therefore also strongly relevant. It is thus difficult 
to treat the effects of such randomness systematically within a controlled 
renormalisation group calculation. 

In this paper, by contrast, we study the effects of quenched randomness 
coupling to the local energy density on systems whose pure versions exhibit 
fluctuation-driven first order transitions. These are transitions which are 
expected to be continuous on the basis of a mean field analysis, but which are 
driven first order by the fluctuation effects. In terms of the renormalisation 
group, they are often characterised by so-called runaway behaviour, that is, 
the renormalisation group trajectories move out of the region in which the 
original perturbative calculation is valid. That, in itself, does not guarantee 
that the system in question undergoes a first order transition, but often 
it is possible to argue that the trajectories then move into a region where 
mean field theory is applicable, and which may then predict a first order 
transition. The Imry-Wortis argument [0, || should, of course, apply equally 
well to systems exhibiting fluctuation-driven first order transitions. However, 
the advantage of studying these from the renormalisation group point of view 
is that it is possible to analyse them within a controlled perturbative scheme 
and to elucidate the nature of the fixed point which governs the continuous 
critical behaviour of the random system. 

A simple example of a two-dimensional system which exhibits a 
fluctuation-driven first order transition is that of N Ising models coupled 
through their local energy densities. Microscopically this may be represented 
in terms of a lattice model with N Ising spins (si(r), . . . , Sjv( r )) a ^ eacn s ^ e 
r of the lattice. The reduced hamiltonian is 

H = ~KJ2I: Si(r) Si (r>) Si (r) Sj (r')^(r) Sj (r'), (1) 

i r,r f i^j r,r' 

where the sums over (r, r') are over nearest neighbour sites. Such a model is 
self-dual on the square lattice, so that the critical coupling K c may be found 
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exactly. In the absence of randomness, the renormalisation group equations 
on the critical surface have the form 



dg/d£=(N-2)g 2 + 0{g 3 ). 



(2) 



For N = 2, this vanishes, as expected since this case corresponds to the 
Ashkin- Teller model which exhibits a line of fixed points labelled by <?.|| 
When N > 2, however, initially small positive values of g flow out of the 
region of validity of the perturbative equation (|2|). When ([J) is analysed 
within mean field theory, the quartic term in the free energy remains positive 
if g is small, indicating a continuous transition, but, for sufficiently large g, it 
changes sign so that the mean field transition becomes first order. Since the 
renormalisation group indicates that, no matter what the initial value of g, it 
should ultimately renormalise into this region, it implies that the transition 
should be first order for all g, and is therefore of a fluctuation-driven nature 
for small g. In fact, on the critical surface, this model when expressed in 
terms of Ising fermions is nothing but the Gross-Neveu modelfl(|], which is 
believed to be massive for N > 2, corresponding to a finite correlation length. 

We now consider adding quenched randomness which couples to the local 
energy density. This may be done in a variety of ways, but, in order to 
focus on the universal properties of such a coupling, let us first rewrite (|l|) 
in a continuum notation in terms of the local energy density E^r) of each 
Ising model. The hamiltonian density, close to the critical point of the pure 
system, may then be written 



where H c is the fixed point hamiltonian and t is the temperature deviation 
from the critical point. Quenched randomness is now added by allowing 
t — > t + St(r), where 5t(r) = and St(r)5t(r') = A5(r — r'). Introduc- 
ing n replicas a = 1, . . . ,n and averaging over a Gaussian distribution, the 
replicated hamiltonian density is 



Note that each term has a well-defined behaviour under the duality operation 
under which E" 1 reverses sign. The self-dual critical point therefore remains 



H = H c + tY,E l -gY<E l E, 



(3) 



(4) 




4 



at t = in this parametrisation. It is possible to consider replica-coupling 
terms which break this duality symmetry, but they are all irrelevant close to 
the pure decoupled fixed point. 

The perturbative renormalisation group equations for the couplings follow 
using standard methods. In general for a perturbed hamiltonian density of 
the form TC = 7i c + J2i9i&i, they have the formfTT| 



dgk/di = y k 9k - J2 c ijk9i9j + 0(g 3 ), (5) 
hi 

where yi is the eigenvalue at the unperturbed fixed point, and is the 
coefficient of $^ in the operator product expansion of $j with $j. In the 
present case, these are very easy to work out. Both the interaction terms in 
(H) have a similar form, and are in fact special cases of a very general model 
of Nn interacting Ising models, with a hamiltonian density 

Nn 

H = H C + ^ G pq E p E q . (6) 

The terms with p = q are excluded since the operator product expansion 
of E p with itself in the Ising model yields only the trivial identity operator. 
Normalising the energy density so that E p ■ E v > = 5 PP >, the required terms in 
the operator product expansion are then 

(E p E q ) ■ (E p iE q i) = 5 pp i(E g E q i) + permutations + • • • , (7) 

from which follow the general renormalisation group equations 

dG pq /d£ = -4(1 - 5 pq ) £ G pr G rq + 0(G 3 ). (8) 

r 

Specialising these to the case at hand, we then find, in the limit n — > 0, 
the flow equations 

dgjdt = 4(N -2)g 2 -8gA + --- (9) 
dA/de = -8A 2 + 8(iV-l)^A + --- (10) 
dtjdt = t(l-AA + A(N -l)g) + ---. (11) 

The last equation in fact follows from the second by the renormalisation group 
version of the Harris criterion [O]. For the case N = 2 these equations are 
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equivalent to those found for the random Ashkin- Teller model by Dotsenko 
and Dotsenko[14|. Remarkably, the flows in the critical surface obtained by 
solving these equations may be found in closed form for general N: 



g = const.(A/^)( JV - 2 )/ iV e- 2A /^. (12) 

These are shown in Fig. [I] for the case N = 3. For g = we find that the 
randomness is marginally irrelevant, consistent with the well-known case of 
the random Ising modelfll^, [U|. In the absence of randomness, the flows for 
g > run away to the first order region, as discussed above. However, for any 
non-zero randomness the trajectories eventually curl around and approach 
the fixed point corresponding to N decoupled pure Ising models. Of course 
equations (pj- |TT|) are strictly valid only inside the perturbative region where 
the initial values of the parameters are small, but it is reasonable to expect 
that the topology of the renormalisation group flows should persist at least 
in some finite region around the origin. This topology has two important 
consequences: first, as dictated by the Imry-Wortis argument, the transition 
has become continuous, and secondly, the asymptotic critical behaviour is 
that of the pure Ising model, similar to those cases discussed earlier. In fact, 
by integrating equation ([TT|) for t it may be shown that the specific heat 
has a singularity of the form v41nln(l/t), just as for the random bond Ising 
model[15|], but with an amplitude A oc A~( N ~ 2 ^ N . Flows of the type shown 
in Fig. [I] are very unusual as they violate the c-theoremfT?||. Of course, that 
this can happen is a consequence of the n — > replica limit. 

It is instructive to extend the above analysis to dimensions d = 2 + e, 
since the Imry-Wortis argument leads to no definite conclusion in that case. 
The perturbative renormalisation group equations become 

dg/d£ = a p g+(A(N -2) + 2b 2 )g 2 - (8 - Ab 2 )gA + ■ ■ ■ (13) 
dAjdt = a p A - (8 - 26 2 )A 2 + 8(iV - l)gA + ■ ■ ■ (14) 
dtjdt = t(d/(2-a p )-(A-2b 2 )A + A(N - l)g) + ■ ■ ■ . (15) 

The linear terms are a consequence of the fact that the specific heat exponent 
a p of the pure model no longer vanishes: this determines the eigenvalues of 
g and A at this fixed point according to the Harris criterion. The parameter 
b is the operator product expansion coefficient appearing in E p ■ E q = 5 pq + 
bS pq E p + • • •. With the energy density normalised in this way, it is universal, 
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depending only on d, but it vanishes for d = 2 as a consequence of duality, 
and is therefore presumably small just above two dimensions. (In d = 4, at 
the Gaussian fixed point, b = 2y/2.) When g = 0, the randomness is now- 
relevant, and the critical behaviour is controlled by a non-trivial random 
fixed point at A = 0(a p ). ]l2j However, when b ^ 0, g is in fact relevant at 
this fixed point, and there exists another, more stable fixed point, which, for 
b 1, is located at g « b 2 a p /lQN, A « a p /8. When g > initially, the 
trajectories move towards larger values of g before eventually curling around 
to finish at this new coupled random fixed point. Therefore this gives an 
example of a fluctuation-driven first order transition in d > 2 dimensions, 
which is converted to a continuous transistion, as in d = 2. However, now the 
critical behaviour is controlled by a new, random, fixed point. Such a fixed 
point, if it persists as far as e = 1, would describe the random Ashkin- Teller 
model in three dimensions, at least for small values of g. 

The above calculation breaks down near four dimensions (if not before), 
due to the proximity of the Gaussian fixed point. As a further example 
of what can happen to a fluctuation- driven first order transition in 4 — e di- 
mensions, consider the well-known problem of the O(N), or N- vector, model, 



with cubic symmetry breaking! 18]. This model has A-component continuous 



spins Si(r), and the replicated hamiltonian density is 

w = *E(^) 2 +«E(-s?)W+ v E(5i) 4 - A E (<s?)W- (is) 

i,a i,a i,j,aj^b 

For u v this may be viewed as a continuous spin version of (|j), with 
g = —u. However, in the absence of randomness this model also possesses 
an 0(n) fixed point (which is absent for d = 2) and a cubic fixed point 
where both u and v are non-zero. The perturbative renormalisation group 
equations may be found from the operator product expansion as above: 

dujdt = eu - 8(N + 8)u 2 + 8NA 2 - 48uv + ■ ■ ■ (17) 
dv/d£ = ev-72v 2 -96uv + --- (18) 
dA/d£ = eA- 16(N -2)A 2 -16(N + 2)uA- 48vA + ■■ ■ . (19) 

When A = these exhibit runaway behaviour to the first order region where 
u and v are large and negative, if the initial value of v is sufficiently negative 
(for N > 4, when the cubic fixed point is in the quadrant with v > 0, this 
requires only that u > and v < initially). However, since A does not 
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enter the flow equation for v, this catastrophe still occurs in the presence 
of randomness. We conclude that quenched randomness does not change 
the order of the transition in this case. This is, of course, quite consistent 
with the Imry-Wortis argument, which does not rule out either behaviour for 
d > 2. 

The case where u < and v > corresponds to the same example of N 
coupled Ising models as before, this time near four dimensions. Solving the 
renormalisation group equations for v near the decoupled Ising fixed point, 
we find (using g rather than u) 

dg/d£ = + 8Ng 2 + 32gA + ■ ■ ■ (20) 
dA/d£ = §eA + 16(iV -2)gA + ---. (21) 

When g = 0, these equations have no perturbative fixed point, despite the 
fact that A is relevant. This is the well-known problem of the random Ising 
model near d = 4, and it is cured in a higher order calculation Jl9|, when a 
term 0(A 3 ) appears on the right hand side of (pi]), giving an 0(e 1//2 ) fixed 
point. However, it may be seen from the structure of the other terms in 
( pO| , pT|) that this cannot cure the runaway behaviour which occurs once the 
coupling g is initially non-zero. We conclude that the fluctuation-driven first 
order transition probably persists in this case. 

As a final example, we may quote the case of the complex O(N) model 
near four dimensions, coupled to a long-range U(l) gauge field, known as 
the Abelian Higgs model for the case N = 1. This was argued long ago 
to undergo a fluctuation-driven first order transition near four dimensions 
for sufficiently small iV.||20|| The effect of quenched random impurities was 
studied by Boyanovsky and Cardyf^], who found that for sufficiently weak 
randomness the first order nature of the transition persists, while for stronger 
randomness the trajectories spiral in towards a new random fixed point, 
corresponding to a continuous transition. 

To summarise, we have given examples of how quenched randomness cou- 
pling to the local energy density converts a fluctuation-driven first order 
transition into a continuous one for d = 2, consistent with the Imry-Wortis 
argument, and of how this may or may not happen when d > 2. 

We conclude with a discussion of the conjecture that all random critical 
behaviour in two dimensions is Ising-like. This is based on numerical results 
on the random Ising modeip^, E3L the Ashkin- Teller model and the 
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4-state Potts model[[7|][] (all of which exhibit continuous transitions in the 
absence of randomness), and the 8-state Potts model (which is first order 
in its pure version.) It is backed up by the interface arguments of Kardar 
et al.||, which suggest that the Widom exponent for the random g-state 
Potts model is independent of q for sufficiently large q. (A similar lack 
of dependence on q has been argued for in the case of random Potts spin 



chains |I3[]; however, their critical behaviour is rather different in nature from 
that of the present case.) 

The conjecture in the case of the Ising model and the Ashkin- Teller model 
(close to the decoupling point) agrees with the results of a perturbative renor- 
malisation group analysis [|14|, and with our analysis above: the renormalisa- 
tion group trajectories curl around and end up at the Ising fixed point, giving 
Ising exponents, but with logarithmic (or log-log) modifications. However, a 
similar analysis J27], applied to the random g-state Potts model for q > 2 
indicates the existence of a new random fixed point whose critical exponents 
are not Ising-like, but depend on q. This analysis is valid only when q — 2 



is small, but is consistent with earlier renormalisation group results |26j for 
g=3. 

We have not been able to resolve this discrepancy, but would venture a 
few remarks which, in fact, may seem to confuse the situation further: 

1. The perturbative renormalisation group arguments work with a dis- 
tribution of randomness which is self- dual , corresponding on the lat- 
tice, for example, to an equal distribution of strong and weak bonds 
of strengths K and K* which are dual to each other. Within the per- 
turbative scheme, this is justified, since it may be argued that weak 
randomness which violates self-duality is irrelevant in the renormalisa- 
tion group sense. 

2. However, the interfacial analysis of Kardar et al.||, which treats hor- 
izontal and vertical bonds on a quite different footing, cannot, by its 
nature, respect the duality properties of the model. Indeed, these au- 
thors find it necessary to include negative bonds in the model to access 



1 There is also claimed experimental evidence for the 4-state Potts model [£5| - however 
the randomness discussed there would appear to favour one sublattice rather than another, 
and therefore should couple to the order parameter, corresponding to the random field 
problem. 
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their zero-temperature fixed point, which are excluded in any self-dual 
formulation of the problem. The only zero-temperature fixed point in 
the self-dual random model is the percolation point. 

This leads to the picture that the critical behaviour controlled by a zero- 
temperature fixed point, discussed by Kardar et al.||, and that found 
in the perturbative renormalisation group of LudwigpTI], are simply 



different and correspond to strong non-self-dual randomness and to 
weak self-dual randomness respectively. However, the numerical results 
for the q = 4 and q = 8 Potts models, which appear to find Ising-like 
exponents independent of q, use self-dual randomness in order to locate 
the critical point precisely. They also consider different strengths of 
randomness, with no appreciable difference in their results. 

4. One would expect critical behaviour controlled by a zero-temperature 
fixed point to exhibit hyperscaling violation, as in the random field 
problem. However, the exponents found in the numerical work for q = 8 
are consistent with hyperscaling, that is, with a conventional, finite- 
temperature fixed point as found in the perturbative renormalisation 
group approach. 

Whatever the resolution of this problem, it cannot be that all the uni- 
versal properties of the random g-state Potts model are independent of q, 
even if the exponents are. This is because this critical point separates a q- 
fold degenerate ordered phase from a nondegenerate disordered phase, and 
this degeneracy must reflect itself in the fluctuation contribution to the free 
energy near the critical point, even if the exponents are Ising-like. This may 
be seen in the example of N coupled Ising models discussed in this paper: 
although the relevant critical fixed point is Ising-like, it in fact corresponds 
to N decoupled Ising models, not just one. This will reflect itself in uni- 
versal amplitude ratios which involve the free energy. However, because of 
the expected logarithmic corrections, these may be difficult to analyse from 
numerical data. A cleaner test should be through the value of the effec- 
tive central charge (which measures the finite-size scaling behaviour of the 
quenched free energy |29|, j30[|). In our example, this is c = ~N, and does 
depend on N. It would be very interesting to compute this for the random 
g-state Potts model. To how many decoupled Ising models does the ran- 
dom g-state Potts model correspond at criticality, if indeed its behaviour is 
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Ising-like? 
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Delta 




Figure 1: RG flows on the critical surface for the case N — 3. Several 
trajectories are shown, for different initial values of the randomness A at 
a fixed value of the coupling g between the Ising models. Although they 
initially flow towards strong coupling, they eventually curl back towards the 
pure Ising fixed point. 
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